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ABSTRACT 

We discuss orientational ordering in monolayers of solid hydrogen in view of recent exper- 
imental findings in NMR studies of (ortho) c -(para)i_ c -hydrogen mixtures on boron nitride 
substrate. Analysis of the temperature-concentration behavior for the observed NMR fre- 
quency splitting is given on the basis of a two- dimension (J = 1) C -(J = 0)i_ c -rotor model 
with the quadrupolar coupling constant To = 0.50 ± 0.03 K and the crystalline field ampli- 
tude Vq = 0.70 ± 0.10 K derived from experiment. The two distinct para-rotational (PR) 
short-range ordered structures are described in terms of the local alignment and orientation 
of the polar principal axis, and are shown to be due to the interplay between the positive 
and negative crystalline fields. It is shown that the observed below the 2D site percolation 
threshold c p = 0.72 local structures are rather different from the ferromagnetic-type PR 
ordering suggested earlier by Harris and Berlinsky. 
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I. INTRODUCTION 



Careful studies at low temperatures of the thick films (from 2 to 12 monolayers) [1,2] and 
monolayers [2,3] of ortho-para hydrogen on boron nitride (BN) substrates revealed at low 
ortho-i?2 concentrations new short-range frozen structures. Besides the analogue of known 
[4-13] quadrupolar glass (QG) phase, that emerges [2] in monolayers below the concentration 
c p = 0.72 , which is apparently close to the site-percolation threshold [14] in honeycomb 
lattice, the uncommon para- rotational (PR) phases denominated by PR-A and PR-B phases 
were discovered [2,3]. They are demarcated by a crossover temperature Tx Cxp \c) at which 
the NMR frequency splitting passes through zero (see 2D diagram in Fig.l). Instead of the 
Pa3 structure known in the bulk hydrogen, the herring bone (HB) and pinwheel (PW) 2D 
long-range ordered structures have been the subject of scientific interest since 1979, when 
Harris and Berlinsky made [15] their famous mean-field theory predictions. Meanwhile 
thorough experimental studies on grafoil [16,17] and BN [2,3] substrates registered only the 
PW orientational order at sufficiently high concentrations, i.e. above the percolation limit 

Cp. 

Analysis is given within the scope of the site-disordered microscopic 2D (J = 1) C (J = 
0)i_ c -rotor model, which is introduced on the basis of the 3D-rotor analog developed earlier 
[18-21] for an indepth study of the QG phase. We will give a microscopic explanation of 
the observed temperature-concentration behavior for the orientational local-order parame- 
ters related to the NMR line shapes. We will show that the PR-A and PR-B short-range 
correlated structures are due to the interplay of the frustrated o-i/2-molecular exchange 
interaction with the molecule-substrate interaction. 

II. MICROSCOPIC DESCRIPTION 

Description of the orientational degrees of freedom of the site-disordered ortho-para- 
hydrogen system with pure electrostatic quadrupole-quadrupole (EQQ) intermolecular in- 
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teractions has been discussed extensively within context of the 3D QG problem [6,8,12,18]. In 
general, the thermodynamic rotational states of a given ortho-molecule located at a site i are 
characterized by the a second rank local tensor that has only five independent components: 

A 2 

the three principal local axes (given by vector Lj) , and the alignment Oi =< (1—3 J zi /2) > T 

A 2 A 2 ^ 

and the eccentricity rji =<J xi — J y i>T defined with respect to Lj axes. (Here J ai stands 
for the angular-operator rotational moment of a given ortho-molecule located [8] at site i, 
and < ... >t refers to a thermodynamic average at temperature T.) A thermodynamic de- 
scription, given in terms of the local molecular fields e a i, and e V i conjugate to the local order 
parameters and extended by the crystalline field hi can be introduced [18,19] on the basis 
of the local-order-parameter fundamental equations, namely 

3cosh(v / 3£ r , l /2T) . 

<7j = 1 = — ; la 

2 cosh( v / 3%/2T) + exp[3(e ai +/ii)/2T] 

= v/3sinh(y/3 £rp /2r) 

Vl 2cosh(\/3^/2T) +exp[3{e ai +hi)/2TY 

These equations follow from the conditions of local equilibrium [10] and they are shown [19] 

to be consistent with the density-matrix representation [7]. In what follows, we restrict our 

consideration to a reduced set of local order parameters {Lj, <jj} with rji = that corresponds 

to the so-called "powder approximation" common [5] in NMR theory applications. This 

description ignores the local field e V i, conjugate to the local eccentricity rji, and Eq. (1) 

(where e a i = E{ and e V i = 0) is therefore reduced to 



1 - Oi 



, , Vt exp[~ ^ £t with £ i = -J2 Ji i°i c r ( 2a ) 

j^i 

Here the effective exchange interaction Jij and the crystalline field hi are given by 

Jij = -^r P 20 (L zi )P 20 (L zj ) and h t = ^V P 20 (L zi ). (2b) 

This reduced mean-field description formally follows from the truncated 2D Hamiltonian 
given for iV quantum rotors with z neighbors placed in the plane, namely 

N z N „ 9 

A \ ^ \ ^ A ^ x ^ A A O ^ 

Hn= ~ 2^2^ 3*3 ° l a i° iC i ~ hi 0i Ci wlth ° i= 2 Jzi ' ^ ' 

i=l j^i i=l 



Here q is a random occupation number whose mean, given by the configuration^ average, 

A 

is the concentration: c =< q > c ; J zi is a ^-projection of the angular momentum operator 
in the local principal coordinate system. In turn, r stands for the EQQ coupling constant 
and Vo is the crystal-field amplitude; P 2 o(-^zi) = (3 cos 2 9j — l)/2 where 0; is the polar angle 
of the principal molecular axis Lj. 

In bulk HCP solid ortho-para-hydrogen the coupling constant T and the amplitude V 
are well established theoretically that is not the case of commensurate a/3 x \/3 solid mono- 
layers. The magnitudes of T = 6Q 2 /(25-Rq) (Q is the molecular electrostatic quadrupole 
moment and i?o is the nearest-neighbor molecular separation) calculated for graphite and 
BN substrates are, 0.534 K [15] and 0.470 K [2], respectively. The approximate estimates 
for the crystal-field amplitude \V (cxp) \ = 0.6 - 0.8 K for grafoil [23] and |V/ (cxp) | « 0.6 K 
for BN [17] were derived from the observed NMR line shapes. As seen from Fig.l, the ex- 
perimentally established position of the order-disorder boundary is consistent for the cases 
of grafoil [16] and BN [2] substrates (shown by crosses and triangles in Fig.l, respectively). 
On the other hand, the PW-PR boundary (restricted by points a and b in the insert of 
Fig.l) exists only for positive crystalline fields. These yield the following fundamental model 
parameters 

T = 0.50 ± 0MK and V = 0.70 ± 0.10K (2d) 

in order to specify our estimations based on the 2D (J = 1) C (J = 0)i_ c -rotor model given 
in Eq.(2). 

III. MACROSCOPIC DESCRIPTION 

The phase diagram for the pure J = 1 rotor system on a 2D triangle solid lattice was 
scaled [15] by Harris and Berlinsky in terms of the EQQ coupling constant r and the crystal- 
field amplitude V of both signs (see insert in Fig.l). The PR phase was postulated by a single 
ferromagnetic-type structure that can be given as {6j = 0, <7j = a }. Moreover, one can see 
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that equation for the PR alignment do obtained by minimization of the relevant free energy 
(see Eq.(17) in Ref. [15]) is equivalent to Eq.(2) with adopted z = 6, = e = —9T a , 
Jij = — 3r /2 and hi = ho = 2Vo/3 that corresponds to a ferrorotational-type (FR-type) 
local-structure given by 0, = Qj = and Oi = er . 

A description for the long-range orientationally disordered, but locally correlated PR-A 
phase is introduced through the short-range order parameter <ta(c, T) =< (Ti(T) >^ RA \ 
where a configurational average is limited by the temperature-concentration PR-A region 
shown in Fig.l. Application of this average procedure to both the sides of Eq.(4) can be 
presented in the following form, namely 



PR-A: = exp 

1+2(74 



Vq _ 3 , Sia + Seia , 9 ,e 2 A^2 
T 2 [ f ' 8 [ ^F } 



(3a) 



Unlike the case of the QG, we assume here that fluctuations of the local alignment (or the 
quadrupolarization) are small. The same is referred to the crystalline field given by the 
mean h 1A = 2V /3. The local fluctuations of the molecular field are introduced through 
the mean E\a = Jiazc'^o a (with J XA = ~ 371^0/2 and with z — 6) and the variance 
e\ A = 3 l 7 2 4 4 2;c(l — c)(l — a A ) 2 /8T 2 , and are estimated [25] within the Gaussian distribution 
justified in Ref. [19]). With account of the Zeeman-field local polarization effects [18] given 
by the mean 5e±a = —J$ a zco ~a{]- — °~a)/T, we have analyzed [24] a concentration behavior 
of the observed [23] NMR frequency splitting given in Fig. 8 in Ref. [2] for T = 0.65 K 
and T = 0.546 K. Analysis is given with the help of Eq.(2) where the polar-principal-axis 
correlation parameters 7^ and 72A, namely 



7 L4 =< P2 (L zi )P 20 (L zj ) >£ RA) and l2A = \J < Pi (L zi ) P| (L zj ) >% RA) (3b) 

are treated as fitting parameters. In the particular cases of the FR-type (Gj = Qj = 0) and 
AFR-type (0; = 0, Qj = n/2) locally correlated structures are characterized by 71 = 7 2 = 1 
and 71 = —72 = —1/2, respectively. For the PR-A phase we have derived [22] 7^ p ' ) = — 1/3 
and 7 2 \4 P ' ) ~ 0.75. This local structure is in a way similar to that in the PW phase modified 
by orientations of in-plane rotors which show out-of-plane orientations. 



As seen from Fig.l, the long-range disordered PR-B phase is stable at low temperatures 
(T < Vo) and low concentrations (c < c p ) where the site-dilution effects are expected to be 
more pronounced than in the PR-A phase. The short-range orientational arrangement results 
from the interplay between the random EQQ coupling and the random negative crystalline 
fields. Adopting for the latter a Gaussian distribution, and taking into account its variance 
h 2 B (with the mean his = 2Vo/3) one finds, after elaboration [25] of the configurational 
average in Eq.(2a), the effective amplitude of the crystalline field can be introduced as 

V(c,T) = V (l - x - U ) for T^T x ^\c). (4) 

Here T x ( exp \c) is a crossover temperature between the PR-A and the PR-B structures that 
provides a reconstruction of the local order from a^ xp ^ (c, T) > to a^ xp ^ (c, T) < (shown by 
dashed line in Fig.l). The explicit form in Eq.(4) follows from /i 2 b — (< A/if >(^'- R - B )) 1 / 2 = 
2(2V T X ) 1 / 2 /3 where T x is approximated by the observed PRA-PRB boundary. An analysis 
of the observed PR-B quadrupolarization o B is given through averaged Eq.(4), namely 

PR-B:ln(i^) + f (1 - - " 4) = 0. (5) 

Treating the PR-B phase as a precursor of the 2D QG phase, we have omitted in Eq.(5) 
all molecular-field local ordering effects. Similar to the QG case, we have therefore adopted 
< Si >^ RB ^= 5eib employed in Eq.(3a) for the PR-A phase. Analysis of the available 
experimental data for c = 0.44 (with T'i cxp " ) = 1.64 K , see Fig. 12 in Ref. [2]) on temperature 
dependence of the short-range orientational order parameter in the PR-B given with the help 
of Eq.(5) results [22] in, approximately, 715 = 0, 72B = 1, that in a way is characteristic for 
the QG local order. 

The observed order parameters <r^ xp - ) and a^ xp ^ vanish at a certain crossover temperature 
T x (c) associated with the PRA-PRB boundary T x exp \c) (shown by the dashed line Fig. 1). 
For concentrations c < c p , this boundary can be therefore deduced from the conditions 
(Ja(c, T x ) = <Jb(c,T x ) = 0. To satisfy the boundary observation conditions, the interplay 
between the fluctuating crystalline and Zeeman-type molecular fields for T <T < T x is 
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made implicit in the form 8V (T - T ) - 9e|(c,T) = 0, where T (= 9hl/8V ) plays a role 
of T x when the competing fluctuations the EQQ field are ignored. The variance h 2 was 
studied [20] in detail for the 3D disordered PR phase in o-p-H 2 systems for temperatures 
0.80 < T < 4.9 K . As seen from the right insert in Fig. 2, unlike the mean of the crystal- 
field hf D \ its variance depends strongly on the overall concentration, i.e. , h 2 ~ c(c M — c) 
and disappears at the highest concentration for the 3D QG state, cm = 0.55 (for the 3D 
phase diagram see Fig. 2 in Ref. [19]). In the 2D case cm is very close to the threshold 
concentration c p . Therefore, we adopt V 2 ~ c(c p — c) that reduces the aforegiven boundary 
observation condition to the following cubic equation 



Ti?-T„(c)T|-('| 72 r ) 4 £^ = with T = m 



Af (1 - f ) 



2 



(6) 



Treating A as an adjustable parameter characterizing a scale of the crystal-field fluctuations, 
we analyze in Fig. 2 the physical solution T x of Eq. (6) by comparing it with the observed 
PRA-PRB boundary. Taking into account the above analysis for PR-B and PR-A phases, we 
adopt 72 = 1 as a typical value. As seen from Fig. 2, the idea that the disordered PRB phase 
is constructed from mostly disordered "in-plane" rotors is corroborated by experimental 
observations [2] . On the other hand, our consideration of the reduced orientational degrees 
of freedom fails to give quantitative descriptions above c = 0.45. The unphysical value 
72 > 1 deduced from experimental data in Fig. 2 signals the existence of ignored local order 
parameters {e.g. q a ^ a 2 ), which (similar to the case of the 3D QG given in Eq.(6) in Ref. 
[18]) can play an appreciable role near the crossover temperature. A complete analysis should 
be given beyond the "powder approximation" and based on fundamental order-parameter 
equation (1) given for both the local alignment and the eccentricity. 



IV. CONCLUSIONS 

We have discussed the short-range orient at ionally correlated structures discovered [2,3] in 
monolayers of (o-H 2 ) c {p-H 2 )i-c on BN substrate. Analysis of the temperature-concentration 
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behavior for the observed NMR line shapes, related to the short-range order parameter 
<r(c, T), is given on the basis of the 2D (J = 1) C -(J = 0)i_ c -rotor model, for which a 3D 
analog was employed earlier [18,19] for the QG problem. In the current study the focus 
is on the nearest-neighbor correlated structures observed by the NMR spectroscopy in the 
orientationally disordered phases (shown in Fig.l). In spite of the fact that fundamental 
order-parameter equations are consistent with the corresponding Eq.(17) in Ref. [15], the ob- 
served PR-A and the PR-B structures are rather " antiferromagnetic" than "ferromagnetic" 
as suggested in Ref. [15] for a unique PR phase. This conclusion follows from our anal- 
ysis of the observed [3,2] macroscopic quadrupolarizations oa(c, T) and <7b(c, T) adjusted 
through the polar-axis correlation parameters given in Eq.(3b). We have shown that the 
short-range correlated PR-A phase is driven by positive crystalline fields, for which ther- 
mal and spatial fluctuations overwhelm those of the short-range EQQ interactions. With 
decreasing temperature, the interplay between the EQQ coupling and the crystalline fields, 
which are both sensitive to the site-dilution and thermal-fluctuation effects, results in the 
PRA-PRB boundary Tx° xp \c), along which both the quadrupolarizations are zero (see anal- 
ysis in Fig. 2). The low-temperature PR phase, denominated as the PR-B phase, is driven 
by negative crystalline field given near the boundary Tx Cxp ^ in Eq.(4). Similar to the case 
of the QG phase (see Fig.4 in Ref. [19]), this phase is expected to be richer than the PR-A 
phase, and more order parameters are therefore needed to give a complete description of the 
observed aJ xp) (c,T). Unfortunately, this data (such as on Aq B =< af(T) - rf(T) >(f i?B )) 
is at present time not available experimentally. 
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Figure Captures 

Fig.l. Phase diagram for site-disordered monolayers of ortho-para-hydrogen (o-if 2 ) c (p- 
H 2 )i- C mixtures. The symbols refer to observed changes in the NMR line shapes reported 
in the literature: crosses for hydrogen monolayers on graphite, Ref. [16]; open triangles 
for commensurate hydrogen monolayers on BN, Ref. [2]. The solid symbols refer to NMR 
studies of Ref. [2]: solid circles, transitions to the quadrupolar glass (QG) state; diamonds, 
transitions to the hindered rotor (HR) state. The inverted triangles refer to the vanishing of 
the small splitting of the NMR lines in the para-rotational state. Insert: theoretical phase 
diagram from Fig. 2 in Ref. [15]; a and b are the tricritical points [15], and c is the minimum 
in the observed PR-PW transition temperatures. 

Fig. 2. Pararotational A-B crossover temperature against concentration. The symbols refer 
to the experimental points represented in Fig. 1: solid inverted triangles, vanishing of NMR 
doublet, Ref. [2]; open triangles, onset of PW state in the anomalous upturn region of the 
phase transition boundary, Ref. [2]. The lines designate solutions of Eq. (6) for an adjustable 
parameter A = 2.3. Other parameters shown include the fitting parameter 72 . Insert: Left; 
Concentration dependence of the effective crystalline field at distinct temperatures derived 
from experiment [2] using Eq.(4); right, the variance of the crystalline field in bulk ortho- 
para-H 2 (squares, from insert in Fig. 5 of Ref. [20]). 
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